Coulomb gauge approach to scalar hadrons 
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Abstract. The Coulomb gauge model, involving an effective QCD Hamiltonian in the Coulomb 
gauge, is applied to scalar hadrons. Mass predictions are presented for both conventional qq meson 
and qqqq tetra-quark states. Mixing matrix elements between these states were also computed and 
diagonalized to provide a reasonable description of the scalar spectrum below 2 GeV. 
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INTRODUCTION 

The scalar hadron spectrum has been a puzzling problem that has attracted wide interest. 
In particular for the isoscalar J PC = ++ channel, below 1 GeV the /o(600) or o has 
had a somewhat confusing history while above 1 GeV the existence of exotic glueballs, 
hybrid mesons and tetra-quark systems has yet to be established. In this paper we 
apply the Coulomb gauge (CG) model to scalar systems to provide further insight 
regarding their structure. In the next few sections we discuss the model, summarize 
selected previous results and present numerical predictions for the scalar spectrum. More 
comprehensive details can be found in Refs. [1, 2, 3, 4, 5, 6, 7, 8]. 

COULOMB GAUGE MODEL 

In the Coulomb or transverse gauge the exact QCD Hamiltonian is 

Hqcd = H q +Hg + H q g+Hc (1) 

H q = J "«/x^(x)[-ia-V + j3inpP(x) (2) 
H g = l -Jdx[^- l n a (x).^Tl a (x)+B a (x).B a (x)} (3) 
H qg = gJdxJ a (x)-A a (x) (4) 

H C = J dxdyf- l p a {x)K ab (x,y)fp b {y) , (5) 



where *F is the quark field with current quark mass m, g is the QCD coupling con- 
stant, A a = (A^Aq) are the gluon fields satisfying the Coulomb gauge condition, 
V • A a = (a = 1, 2, ...8), IF = -W tr are the conjugate momenta and B" = VxA fl + 
\gf abc K b x A c , are the non-abelian chromodynamic fields. The color densities, p a (x) = 
^(x)T ax ¥(x) +f abc A b (x) -n c (x), and quark currents, J a = ^ (x)aT ax ¥(x), contain 
the standard SU (3) color matrices, T a = 4^, and structure constants, f ahc . The Faddeev- 
Popov determinant, ^ = det(^#), is a measure of the gauge manifold curvature and 
involves the color matrix ^# = V • D with co variant derivative D ab = 8 ab V — gf abc A c . 
The kernel in Eq. (5) is given by K ab (x,y) = (x,a\^^V 2 ^- l \y,b) . The Coulomb 
gauge Hamiltonian preserves rotational invariance, is renormalizable, permits resolu- 
tion of the Gribov problem, avoids spurious retardation corrections, aids identification 
of dominant, low energy potentials and introduces only physical degrees of freedom (no 
ghosts). 

Our model entails two approximations, replace the exact Coulomb kernel with a 
calculable confining potential and use the lowest order, unit value for the the Faddeev- 
Popov determinant, giving the CG model Hamiltonian, Hqg = H q + H CG + H qg + H GG 

Hf G = I J dx[n fl (x)-n a (x)+B fl (x)-B a (x)] (6) 

H£ G = - X -jdxdyp a (x)V(\x-y\)p a (y). (7) 

A Cornell type potential, V(r) = —a s /r+ or, is used for the confining kernel with 
previously determined string tension, o = 0.135 GeV 2 , and a s = 0.4. 

Next, hadron states are expressed as dressed quark (anti-quark) Fock operators, fij^ 

(D^), with helicity, A = ±1, and color = 1,2,3 acting on the Bardeen-Cooper- 
Schrieffer (BCS) model vacuum, \C1) (see Refs. [3, 7] for full details). The qq meson 
state is 

^ JPC) = /(^<£(kK^(k)Z)^(-k)|a). (8) 
For the tetra-quark system the wave function ansatz 

im/PCx _ f gjU d( l2 dq 3 JPC 

} -7(2^(2^(2^ 11 W4(qi ' q2,q3) ( ) 

is adopted with quark (anti-quark) cm momenta qi, q3 (q2, q4). The expression for 
the matrix R%£g 4 depends on the specific color scheme selected [6, 8]. For the color 
singlet-singlet scheme, [(3 <g> 3)i <8> (3 <8> 3)i]i, where the qq pairs, A and B, couple 
to color singlets, = S&^Stf^- This gives the lowest mass among the four 

color representations. The spin wave function part is, (^Ai^sa^-a) (jj^3^4\sb^b) 
(saSb^a^b\J^a + ^b), a product of Clebsch-Gordan coefficients where J is the total an- 
gular momentum, Sa = Si + S2 and Sb = S3 + S4. For all scalar tetra-quarks states the 



orbital angular momenta are zero, consistent with the lowest energy state. A Gaussian 
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radial wavefunction is used (see [8] for details) f(qA, <3b, qi) = e "a "b a i 7 with vari- 
ational parameters OCa = OCb and Gt] determined by minimizing the tetra-quark mass 

Mjpc = {^ JPC \H CG \^ JPC )=M self + M qq +M qq +M qq +M annih , (10) 

which was previously calculated [6, 8]. The respective contributions are the q and q 
self-energy, the qq, qq and qq scattering, and the qq annihilation. 

APPLICATION TO SCALAR HADRONS 

The most recent particle tabulation [9] lists nine isoscalar ++ states: /o(600),/o(980), 
/o(1370),/o(1500),/o(1710),/o(2020),/o(2100),/ (2200) and / (2330), of which the 
last four are not included in the summary table. Our model has previously predicted that 
there is one gg glueball around 1700 MeV [4] and two hybrids: a ring at 2135 MeV, 
where nn = -^(uii + dd), and a ssg at 2140 MeV. In the pure quark sector, we also 

calculated [3] two conventional qq states which were orbital p-waves: a nn at 848 MeV 
and a ss at 1297 MeV. Most recently we have predicted [8] several tetra-quark states: two 
nnnn at 1282 and 1418 MeV and two nnss at 1582 and 1718 MeV. There will also be a 
ssss as well as radially excited nn and ss states near to above 2 GeV which we have not 
calculated. There are three model corrections which should be included and will modify 
this predicted spectrum. The first is chiral symmetry which will significantly lower one 
of our predicted nnnn states corresponding to the kk quantum number channel (this is 
the 1418 MeV state). This effect is discussed further below and will be incorporated in a 
future analysis. The second is to perform a resonance or scattering calculation to obtain 
the imaginary part of the pole position or width of the state. This will also slightly affect 
the real part of the pole position, or resonance mass. We have begun such an analysis [10] 
and will report subsequent developments elsewhere. The third effect is mixing which we 
now address. 

We only treat qq meson and tetra-quark mixing since in calculating quark-hybrid and 
quark-glueball mixing matrix elements with our model Hamiltonian, the former are per- 
turbative, and thus expected weak, while the latter entirely vanish (mixing must proceed 
via higher order intermediate states). This suggests that glueball widths might not be 
large, as typically expected, consistent with a recent theoretical prediction [10]. Mixing 
with gluonic states clearly merits further study which we plan to address in a future 
analysis. Using the notation, \qq > and \qqqq > for \m JPC >, the mixed state is given by 
\J PC ) = a\nn) +b\ss) + Ci\nnnn)i + di\nnss)i for i = 1,2. The coefficients a,b,Ci and d[ 
are determined by diagonalizing the Hamiltonian matrix, in which the meson-tetra-quark 
off-diagonal mixing element is (only H^ G contributes) M = (qq\H^ G \qqqq) , where \qq) 
is \nn) or \ss), and \qqqq) is \nh~nn) or \nnss). There are 12 off-diagonal matrix elements 
however three, (ss\H GG \nn) and (ss\H^ G \nnnn)i, vanish and four, i(nnnn\H GG \nnss}i, 
are numerically very small. The remaining mixing matrix elements are, (nn\H GG \nnnn) 
(nn\H GG \nnss)i and (ss\H GG \nnss)i. For our model Hamiltonian, there are two types of 
mixing diagrams illustrated in Fig. 1. Because of color factors, nonzero mixing only 



FIGURE 1. Diagrams for the meson, tetra-quark mixing term. 



exists for qq annihilation between different singlet qq clusters. The first diagram gives 
Mi = \J Jq 1 Jq 2 Jq3y(fc)^(q 1 )^(-q4)F lll4 (q 1 ,q 2 ,q3)4»i^ 4 (-2q 4 ) , (11) 

with q 4 = -qi - k, k = q 2 + q 3 , F Xl x A (qi , q 2 , q 3 ) = ^ (q^ (q^^J^ fai > ^ q 3 ) 
and dressed, BCS spinors ^ and Yx . The effective confining potential in momentum 
space is V(k). The second diagram yields 

M 2 = rfqi^q3V(fc)^(q4)^(-q 1 )F AlA4 (q 1 ,q 2 ,q 3 )^(2q 1 ) . (12) 

The two Hamiltonian parameters in our model were independently determined while 
the wavefunction parameters were obtained variationally. Because we seek new model 
masses, the unmixed variational basis states need not be ones producing a minimal, 
unmixed mass, so we selected one of the variational parameters, a/, to provide an 
optimal mixing prediction and then studied the mixing sensitivity to this parameter. 

For the ++ tetra-quark state, the spin of the two qq clusters must either be sa= sb = 
or sa = sb = 1, and for each the three mixing matrix elements versus «/ are shown in 
Fig. 2. The mixing term is zero when «/ is zero and then increases with increasing «/. 
Note that mixing with ss states is stronger than with nn states. 




FIGURE 2. Mixing matrix elements versus aj. Solid (dashed) lines are for qq spin (1). 



TABLE 1. Mixing coefficients and masses in MeV for isoscalar ++ states. 



nn> \ss> \nnnn>\ \nnnn >2 \nnss >\ \nnss >2 



no mixing 


848 


1297 


1282 


1418 


1582 


1718 


mixing 


783 


1026 


1295 


1466 


1611 


1962 


exp. 


/o(600) 


/o(980) 


/o(1370) 


/o(1500) 


/o(1710) 


/o(2020) 




400 - 1200 


980±10 


1200 - 1500 


1507 ±5 


1718±2 


1992 ±16 


coeff. 


a 


b 


Cl 




d\ 


d2 


/o(600) 


0.947 


-0.070 


0.138 


0.277 


0.006 


0.030 


/o(980) 


0.052 


0.838 


-0.012 


-0.021 


0.122 


0.529 


/o(1370) 


-0.181 


0.026 


0.973 


0.143 


0.002 


0.004 


/o(1500) 


-0.248 


0.030 


-0.187 


0.950 


0.005 


0.009 


/o(1710) 


-0.036 


-0.315 


0.000 


-0.007 


0.901 


0.295 


/o(2020) 


-0.053 


-0.439 


0.000 


-0.006 


-0.416 


0.795 



Using the calculated matrix elements and previously predicted unmixed meson and 
tetra-quark masses [3, 6, 8], the complete Hamiltonian matrix was diagonalized to 
obtain the expansion coefficients and masses for the new eigenstates. For aj = 0.2, 
the results for ++ states are compared in Table 1 to the observed [9] lowest six ++ 
states. Noteworthy, after mixing, the nn meson mass is shifted from 848 MeV to 783 
MeV and the strange scalar meson mass also decreases from 1297 MeV to 1026 MeV, 
now close to the experimental value of 980 MeV Mixing clearly improves the model 
predictions as the masses of the other /q states are also in better agreement with data. 
Figure 3 illustrates the over all improved description that mixing provides for the fo 
spectrum. New structure insight has also been obtained from the coefficients, with the 
predictions that the c//o(600) is predominantly a mixture of nn and nnnn states while 
the /o(980) consists mainly of ss and nhss states. There is general consensus that 
the <t//o(600) state is a %% resonance (pole in the %% scattering amplitude) and thus 
has a predominantly tetra-quark nature. Because the current numerical treatment of our 
CG model neglects chiral symmetry, our mass predictions for tetra-quark states which 
couple to kk quantum numbers will generally be too heavy and should be regarded 
as preliminary. A chiral symmetry preserving RPA variational mixing study will be 
conducted in the future which should yield a lighter scalar tetra-quark mass since our 
approach is formally equivalent to the Schwinger-Dyson treatment which previously 
predicted the o is a nn resonance [11]. 

SUMMARY AND CONCLUSIONS 

Using the established CG model we have calculated qq and qqqq mixing for the low- 
lying ++ spectrum. We find mixing effects are important and necessary for an improved 
hadronic description. Future work will address mixing applications to glueball and 
hybrid mesons to facilitate establishing their existence. 
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FIGURE 3. Predicted unmixed, mixed and experimental fy spectrum. 
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